The "devil's staircase" type phase transition in NaV205 under high pressure 
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The "devil's staircase" type phase transition in the quarter- filled spin-ladder compound NaV2 05 
has been discovered at low temperature and high pressure by synchrotron radiation x-ray diffraction. 
A large number of transitions are found to successively take place among higher-order commensurate 
phases with 2ax2bxzc type superstructures. The observed temperature and pressure dependence of 
modulation wave number qc, defined by 1/ z, is well reproduced by the Axial Next Nearest Neighbor 
Ising (ANNNI) model. The qc is suggested to reflect atomic displacements coupled with charge 
ordering in this system. The experimental fact implies that two competitive inter-layer interactions 
between the Ising spins, i.e., the nearest neighbor Ji > (ferro) and the next nearest neighbor 
J2 < (antiferro) along the c-axis, are intrinsic in this compound. A microscopic origin of the 
inter-layer interaction is not yet known; however, the nearest neighbor interaction Ji > between 
the V2O5 layers is especially interesting. It is very surprising that the phase transition of such a 
complicated charge-lattice-spin coupled system NaV205 can be described by the simple ANNNI 
model. 
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Since discovery of apin-Peierls-like behavior in NaV2 05 
by Isobe and UedaEI, a wide variety of experimental 
and theoretical studies-|]aave beea carried out. The re- 
cent structure analysi^'H, NMRo and convergent-beam 
electron diffractionU experiments have confirmed that 
NaV205 has an orthorhombic unit cell (P-phase called 
hereafter) with the space group of centrosymmetric 
Dg^-Pmmn and lattice constants a = 11.325 A, 6 = 
3.611 A, c = 4.806 A under ambient conditions. Na 
atoms located between V205-layers composed of two- 
dimensionally connected VOs-pKramids in the a6-plane 
play a role of an electron donomfl. All crystallographi- 
cally equivalent V atoms are dressed a half of 3d^ electron 
|'Y4.5+^ so that spins (S = 1/2) form a quarter-filled lad- 
der structure along the 6-axis. At Tc — 35 K, the charge- 
orderin^fl {Y'^-^t,-^'^-^^ ^Y'^^-Y^^), spin-gap formation 
(A = 9.8 meV)ly'Q and lattice dimerization character- 
ized by the modulation wave vector q=(l/2, 1/2, 1/4|) 
(2a X 26 X 4c in real space; Ci/4-phase called hereafter)Ll, 
take place cooperatively. Regarding the charge order- 
ing pattern below Tc, Seo & Fukuyamaa proposed a 2- 



dimensional zigzag type charge ordering model in the ah- 
plane based on the Hartree approximation by taking both 
on-site and inter-site Coulomb interactions into account. 
Recently, Nakao et alU carried out anomalous x-ray scat- 
tering experiments to selectively identify V''+ and 
ions and directly evidenced for the zigzag type charge or- 
der realized in the a5-plane as theoretically proposed. A 
full structure analysis of both atomic displacement and 
charge order modulation with q=(l/2, 1/2, 1/4) in the 
low-temperature^nhase is very recently accomplished by 
x-ray scatteringj'Q. However, any microscopic mecha- 
nism to stabilize the quadrupling of its unit cell along the 
c-axis is not yet known. The pressure application will be 
effective for providing a key information to clarify such 
an origin of the long period lattice and charge modula- 
tion. This is one of purposes of the present high-pressure 
x-ray diffraction experimeate. 

Recently, Ohwada et aZ.t3 have discovered a new high 
pressure phase which has a 2a x 2& x Ic superlattice 
(Co-phase) led by the Ci/4-to-Co phase transition tak- 
ing place at 0.92 GPa (T = 8 K). On the other hand. 
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dielectric measurements precisely performed along the c- 
axistil indicated the existence of an intermediate phase in 
a pressure range of 0.5 GPa < P < 1.0 GPa at tempera- 
ture T = 20 K - 30 K. The second purpose of the present 
synchrotron x-ray diffraction experiments at low temper- 
atures and high pressure (LT-HP) is to clarify such an 
intermediate phase from a microscopic structural points 
of view. 

For an overall survey of superlattice reflection in re- 
ciprocal space, an x-ray oscillation photographic method 
was first employed on Micro Powder Diffractometer 
(MPD)E3 equipped with an imaging plate (IP) and with 
an x-ray CCD systemllj at Photon Factory BL-IB. Wave- 
length was tuned to 0.6888 A (18 keV) with a Si(lll) 
double-crystal monochromator. For high resolution mea- 
surement to obtain more detailed information, a detec- 
tor method was employed on a four-circle diffractometer 
at BL-4C (PF) and BL02B1 (SPring-8). X-ray wave- 
length monochromatized in the same way as that of BL- 
IB (PF) was 0.6888 A at PF and 0.413 A (30 keV) at 
SPring-8. For the LT-HP experiments, a He-gas driven 
diamond anvil cell (DAC) was mounted on a closed-cycle 
He-gas refrigerator. Pressure was generated in a DAC 
using a 4:1 mixture of methanohethanol pressure trans- 
mitting jHcdia and was calibrated from a lattice constant 
of NaClEj enclosed with the specimen in the DAC. The 
sample temperature was monitored with a Au-0.07%Fe 
chromel thermocouple at BL-IB (PF) and a Si-diode sen- 
sor at BL-4C (PF) and BL02B1 (SPring-8) directly at- 
tached to the DAC surface. We used as-grown high qual- 
ity single crystals of NaV2 05 with a typical size of 200 
/iTOx400 fimx50 fim (oxfoxc) grown by the same method 
as previously reportedt3. 
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FIG. 1. Temperature dependence of a position of super- 
lattice reflections (15/2, 1/2 I) in the region recorded on IP 
oscillation photographs at 0.88 GPa. Several intermediate 
phases with I = 1/4, 1/5, 3/16 and 2/11 are systematically 
observed. 

Figure ^ shows the oscillation photographs represent- 
ing temperature dependence of positions of superlattice 
reflections along the direction of [15/2, 1/2, /] in selected 
range of reciprocal space at a fixed pressure of 0.88 GPa. 
All recorded reflections can be completely indexed as in- 
dicated in the figure with the aid of data processing sys- 



tem DENZOEj. With increasing temperature the super- 
lattice reflection (15/2, 1/2, 3/4) disappears at 19 K while 
a new reflection (15/2,1/2,4/5) appears. Furthermore, 
another new reflection (15/2, 1/2, 13/16) appears at 20.0 
K and gradually shifts to the position (15/2,1/2,9/11) 
between 21.5 K and 23.0 K. One can notice that all phases 
have the 2a x 26 x zc type superlattices whose modula- 
tion wave vector along the c-axis can be expressed as qc 
( = 1/z) = 1/4 1/5 3/16 ~ 2/11 in a reciprocal 
lattice units (r.l.u.). The superlattice with qc = 1/6 is 
also observed in a limited temperature range. 
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FIG. 2. Temperature dependence of the superlattice reflec- 
tion experimentally scanned along the [13/2,3/2,/] direction 
with high resolution at (a) 0.92 GPa and (b) 1.04 GPa. (c) 
Peak intensities of both the Ci/4-and Co-phases are shown to 
estimate the pressure distribution. 

Figure |2|(a) displays temperature dependence of high- 
resolution diffraction proflles observed along [13/2, 3/2, 
direction at 0.92 GPa with the counter method. In this 
flgure a few more reflections are newly observed and one 
can clearly see a series of superlattice reflections with 
qc = 1/4, 1/5, 1/6, 2/11 and 3/17 which systematically 
appear and disappear as a function of temperature. Fig- 
ure ||(b) shows the similar data observed at 1.04 GPa. At 
the lowest temperature 11 K, three phases with qc = 1/4, 
1/6 and coexist. As temperature increases, two phases 
with qc = 1/4 and disappear while the other phase 
with qc — 1/5 emerges in a limited temperature range. 
At higher temperatures, the phases with qc = 2/13 and 
1 /7 show up and the three peaks shift to the more compli- 
cated higher order commensurate position. Some phases 
are found to coexist at a certain temperature range as 
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seen in Fig. ||(a) and (b). Such a coexistence may re- 
sult from the pressure distribution in the DAC caused by 
the sohdified pressure media. The pressure distribution 
can be estimated as follows: Fig. ^c) shows pressure de- 
pendence of peak intensities of the C1/4- and Cg-phases 
observed at the temperature where the superlattice in- 
tensity is saturated. The two phases have a clear phase 
boundary at 0.95 GPa nearly parallel to the temperature 
axis. From the pressure range where both phases have 
finite peak intensity at about 1.0 GPa, the pressure dis- 
tribution is estimated as AP = 0.1 GPa. Considering 
thus estimated pressure (AP = 0.1 GPa) and tempera- 
ture (AT = 0.2 K) distributions (defined as the P-T res- 
olution), we obtain the P-T phase diagram of NaV205 as 
shown in Fig. ||. All phases are expressed with a charac- 
teristic modulation wave number along the c-axis qc as 
Cq^, i.e., the C1/4 for qc ~ 1/4 and the Cq for qc ~ 0. 
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FIG. 3. Experimentally observed P-T phase diagram of 
NaV2 05. The hutched area shows commensurate phases un- 
ambiguously identified while the shaded area indicates more 
complicated higher-order commensurate or incommensurate 
phases unresolved by the present P-T resolution, which is rep- 
resented with an ellipse at the left bottom corner. The phase 
boundary experimentally not clearly confirmed is drawn with 
dotted lines. 

We found some more phases whose boundaries could 
not be determined. Such unidentified phase regions are 
represented by dotted lines. The hutched area shows 
the commensurate phases unambiguously identified. The 
shaded area shows more complicated higher-order com- 
mensurate phases or incommensurate phases which we 
could not identify within the present P-T resolution. The 



intermediate phase is no longer observed below P = 0.4 
GPa and above 1.25 GPa (see the inset of Fig. |). At 2.0 
GPa the P-to-Co phase transition temperature decreases 
down to 10 K. 

The observed qc sequences are well understood as the 
"devil's staircase" type sequences theoretically obtained 
from the-simple ANNNI model proposed by Bak and von 
BoehrrO. They developed theoretical calculations with 
two competitive inter-layer interactions between Ising 
spins, i.e., the nearest neighbor interaction Ji > (ferro) 
and the next nearest neighbor interaction J2 < (anti- 
ferro). Within a mean field approximation, they obtained 
a global phase diagram and found extremely complicated 
phase boundaries separating various kinds of commensu- 
rate phases. It is essentially important that there are only 
two stable ground states as T ^ 0, i.e., ttt • ■ ' (ferro, 
all up configuration, q = 0) state for —J2/J1 < 1/2 and 
TTiiTTii-- - (up- up-down-down configuration, q = 1/4) 
state for — J2/J1 > 1/2. Such a transition between the 
two ground states as a function of —J2/J1 was experi- 
mentally observed in NaV205 which undergoes the G1/4- 
to-Co phase transition at 0.92 GPall3. The other charac- 
teristic theoretically proposed for devil's staircase at fi- 
nite temperature is a systematic existence of modulated 
phases with the wave numbers q = 1/5, 1/6, 1/7, 2/11, 
3/17, 3/16, 2/13,... , which is also experimentally ob- 
served in a wide range of temperature and pressure in 
NaV205 as displayed in Fig. ||. Thus the P-T phase dia- 
gram observed in NaV2 05 extremely well resembles the 
theoretically obtained "devil's staircase" with respect to 
the modulation wave number along the c-axis. As for the 
structure in the a6-plane, the 2a x 26 superstructure is 
well maintained. The present experiment directly probes 
the atomic displacements represented by q = (1/2, 1/2, 
qc) but does not directly the charge orderingj. However, 
the previous x-ray scattering experiments Efcl evidenced 
that in the Ci/4-phase the charge order and atomic dis- 
placements are modulated with the same wave vector q = 
(1/2, 1/2, 1/4). In the present phase diagram, therefore, 
it is suggested that the charge order is also modulated 
with the same wave vector in each phase. Such a 2a x 26 
structure is stabilized by the inter-site Coulomb interac- 
tionsO. 

By analogy with the ANNNI model, the devil's stair- 
case type behavior along the c-axis implies any com- 
petitive inter-layer interactions, i.e., phenomenologically 
Ji > and J2 < in NaV205. The inter-layer inter- 
actions must relate to the lattice constant c which has 
compressibility as large as Ac/cq ~ - 0.02 (GPa~^)E3. 
As shown in Fig. ^ (a), Ising spins can be defined by 
the charge disproportionation on the V-O-V rung as S 
= 1 (T) for the V4+-0-V5+ state and S ^ -1 (i) for 
V^"''-0-V*+ state. Local charge arrangement of NaV205 
and Ising spin arrangement of ANNNI model are rep- 
resented in Fig. ^ (b), respectively. It is considered that 
d-electrons in NaV205 are correlated in a long range with 
the Coulomb repulsion. J2 < will be understood eas- 
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ily by the electronic interaction; however, the Ji > 
requires more comphcated mechanisms to overcome the 
Coulomb repulsive energy. The lattice distortion gen- 
erating the polaron or the magnetic interaction between 
the (magnetic) spins can be one of possible origins for the 
attractive interaction between electrons. It is interesting 
to further study a microscopic mechanism of these com- 
petitive interactions and interplays between the lattice, 
charge and (magnetic) spin in this system. 

(a) Ising variable 

^ = i V4-0-V5^ 

(b) 
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FIG. 4. (a) Ising spins correspond to the charge dispropor- 
tionation on V-O-V rung defined in the present system, (b) 
Local charge arrangement of NaV2 05 system along the c-axis 
and Ising spin arrangement of ANNNI system. 
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